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ABSTRACT 

New photometric observations of the hierarchical echpsing TY CrA system were taken 
in the optical with VYS0S6 and in the near-IR with SOFI and REMIR. They are 
the first observations showing the deep eclipse minimum of the pre-main sequence 
secondary in the near-IR. For the first time, the secondary minimum can be reliably 
used in the calculation of the 0-C diagram of TYCrA. By now, the 0-C diagram can 
be studied on a time basis of about two decades. We confirm, that the 0-C diagram 
cannot be explained by the spectroscopic tertiary. For the first time, the light curve 
of the inner eclipsing binary is analysed in both optical and near-IR bands simultane- 
ously. In combination with already published spectroscopic elements, precise absolute 
dimensions and masses of the primary and the secondary component are obtained 
using the ROCHE code. The inclusion of the near-IR data puts strong constraints on 
the third light which is composed of the reflection nebula, the spectroscopic tertiary 
and a visual fourth component. The absolute parameters of the inner eclipsing binary 
agree very well with previous work except of the primary radius (1.46 ± 0.15 Rq) 
and luminosity (40 ± 10 Lq) which are clearly smaller. While the parameters of the 
secondary are well understood when assuming an age of about 3-5Myrs, the primary 
seems considerably undersized. Low metallicity cannot explain the parameters of the 
primary. 

Key words: binaries: eclipsing - stars: pre-main sequence - stars: evolution - stars: 
fundamental parameters. 



1 INTRODUCTION 

The star-forming region R Coronae Australis (R CrA) 
harbors dozens of young intermediate- to low-mass stars 
(and brown dwarfs) of B8 to M8.5 spectral types. They 
have an age of between one and a few million years and 
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a distance of approximate ly 130 pc (for a review see 
iNeuhauser fc Forbrichll2008t ). 

TY CrA is a hierarchical triple, maybe even quadru- 
ple system, embedded in a reflection nebula in the R 
Coronae Australis star forming region. The immediate 
environment of TY CrA h a s bee n studied recently i n 
more detail by ICeers et all (l2007fl. lJuvela et all ll2008h . 



ICurrie fc Sicilia-Aguilar' (|201linSandeU at alj l|201lh and 
Kumar et al. (2011). Two components form amassive eclips- 
ing double- lined spectroscopic binary (eSB2) with an orbital 
period of almost 3 days. While the primary, a B8 star with a 
mass of 3.16 Mq, has already reached the main sequence, i.e. 
initiated stable hydrogen burning, the secondary still resides 
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in the pr e-main sequence ( PMS) stage and has a mass of 
1.64 Mq l|Casev et al.lll99a '). A third spectrosc opic compo- 
nent is in a wide orbit around the eclipsing pair (jCasev et alJ 
ll995l : [ci3rporon et al.|[l99(j) . A vis ual fourth component was 
detected bv lChauvin et~aLl (120031 ) . 

The analysis of eSB2 allows us to derive the most fun- 
damental parameters, mass and radius, as well as to obtain 
effective temperature and luminosity. Thus, the TY CrA 
eSB2 is one of the key objects that can provide precise con- 
straints on the evolution of PMS stars. 

The derivation of effective temperature is a particularly 
important issue. The analysis of the eclipse light curve (here- 
after LC) yields the ratio of the effective temperatures of 
the two components but does not enable to determine in- 
dividual temperatures. Temperatures for both components 
can be derived if a reasonable estimate is used for one of 
them. This step is a major problem in the analysis of PMS 
stars sinc e their effective temperatur e s are not well-known 
a priori (iHillenbrand fc White! 12004 lAmmler et al. I l2005l : 
ISchoning fc Ammlei 20081 ). Because of a significant inter- 
stellar and circumstellar absorption to wards TY CrA, ob - 
served colours have to be used with care JCasev et al.lll99^ . 

TY CrA is one of o nly two known objects (see 
iHiUenbrand fc White! [200i ) which allow us to circumvent 
this problem: the primary has already reached the main se- 
quence, so that its effective temperature can be estimated 
with a reasonable accuracy using the mass - temperature re- 
lation of main sequence stars. The secondary is still located 
in a PMS state. 

The eclipsing nature o f TY CrA was reliably estab- 
lished by !Kardopolov et al] (!l98l! ) using BV photoelectric 
photometry. The authors observed ~0.4 mag deep minima 
consistent with an orbital period of 2.888777 days. With- 
out spectroscopy it was, however, impossible to exclude the 
double orbital period for the system . The first spectroscopi c 
orbit for TY CrA was obtained bv lLagrange et aU ([1993) • 
The authors showed that the radial-velocity (hereafter RV) 
variations of all narrow spectral lines are periodic and con- 
sis tent with the 2.88 8 777 d ays period previously reported 
bv !Kardopolov et al.! (|l98ll ). The width of lines, vsmi < 
6 km s~^, indicated subsynchronous rotation. The authors 
noticed the second system of much wider li nes anticorrelat- 
ing in RV with respect to the narrow lines. ICorporon et al.l 
(|l994! ) confirmed that the broad lines correspond to the sec- 
ondar y component and d erived orbits for both components. 
Later, ICasev et ahl (|l995h found that TY CrA is a SB3 sys- 
tem. The RVs of the tertiary were found to correlate with 
the barycentric RV of the eclipsing pair confirming the grav- 
itational bond. The authors were able to derive a precise 
solution of the circular orbit of the eclipsing binary after 
applying appropriate offsets to the primary and secondary 
data at different epochs. The measurements were insuffi- 
cie nt to derive a unique solution of the tertiary orbit (alse 
see ICorporon et al.! 1 19961 ) which most certainly has a long 
peri od compared to the orbit of the eclipsing binary. 

ICasev et al.l (| 19981 ) obtai ned uvhy photoe lectric pho- 
tometry of TY CrA (see also IVaz et ail 1 1998! ). The mea- 
surements used HD 176423, HD 176497 as comparisons. The 
authors confirmed the orbital period of the eclipsing binary 
of 2.89 days. The primary with spectral type B8 has a mass 
of 3.16 ±0.(32 Ma the secondary 1.64 ± 0.01 Mq. Following 
ICasev et al.! (|l998l ). the age of the eclipsing binary is of the 



order of a few million years according to evolutionary mod- 
els. The LC analysis was challenging mainly because of the 
very shallow secondary minimum in the optical bands with 
a depth of only ~ 0.02 mag. The authors claim that the 
accuracy of the LC is not limited by observational errors, 
which are below 0.01 mag but by ou t-of-eclipse variability 
at a level of 0.05 mag (|Vaz et al.lll998! ). The analysis is com- 
plicated by an offset of the 1992 data relative to the 1993 
data around phase 0.6 which indica tes another type of vari- 



ation, possibly due to stellar spots. ICasev et al.f (|l998h did 



not find a unique solution and pointed out the need for more 
photometric data. In addition, there are effects due to the 
spectroscopic tertiary (light-time effect), the third light de- 
tected in the LC analysis, the asymmetry of the primary 
minimum and the reflecti on nebula to be considered. 

Observations of IChauvin et al.! (|2003| ) with 
VLT/NACO showed a visual companion to TY CrA 
separated by 0'.'294. If the additional component was a 
physical member of the system its spectral type would be 
M4. At such small separation the visual component con- 
tributes to all ground-based photometry and spectroscopy. 
Its light contribution in NIR is about 4% in the J, and 
7.6% in the K passband, respectively. If physically bound, 
the fourth component would have significant consequences 
for the dynamical evolution of the system. 

The aim of this work is to improve the knowledge 
on TY CrA and to solve open issues identified previ- 
ously. Recent work shows that uncertainties can be sub- 
stantiall y reduced when invo lving observations in the near- 
infrared. ICovino et al.! (|2004l ) included near-infrared photo- 
metric data in the analysis of the young eclipsing binary 
RX J0529. 4+0041 and were able to reduce the error bars 
on mass and radius of the components by ~ 80 % compared 
to previous analyses in visual bands. The secondary mini- 
mum was deeper in the infrared band than in the optical 
thus allowing for higher accuracy. Also the determination 
of effective temperature from stellar colours benefits from 
the use of near-infrared o bservations, which are less a ffected 
by interstellar absorption. ISchoning fc AmmleJ (|2008l ) found 
that calibrations based on an optical pass band and an in- 
frared pass band (e.g., [V — K)) give effective temperature 
for young stars which agree best with spectral types. 

The present work for the first time combines NIR light- 
curves (LCs) of TY CrA with observations in visual ba nds. 
The o bservational limitations encountered by Vaz et al.l 
(Il998!) in the visual are overcome by using the advantages 
of wide-field imaging where lots of comparison stars are ex- 
posed simultaneously together with the target. The goal of 
the present paper is to conclusively determine the effec- 
tive temperature of the pre-main sequence secondary and 
to solve uncertainties in previous LCs solutions. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Differential JH photometry 

First JHKs observations for this project were taken during 
two nights (July 30 and September 7) in 2006 with the SOFI 
instrument fsee iMoorwood et al.!ll99^ . and Table [2]) at the 
3.5m ESQ NTT at La SiUa, Chile. The observations were 
scheduled to cover one secondary and one primary minimum 
of the system. 
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Table 1. Parameters of the variable a nd co mparison stars used. The proper motions and {B — V) = 0.8 5(-B — V)t colo rs were taken 
from the Tycho-2 Catalog |Hg g ct al. 20001). The d istance to the system TY CrA was taken from Cliauv in et al.l ||2003| ) and infrared 
colors, (J - K), from 2MASS ijSkrutskie et al.ll2006h : cmp5 = 2MASS 19013729-3649210 = USNOBl 0531.0836104. 



Comp. 


CSC 


HD 


SAO 


lia [mas.y ^] 


IJ.S [mas.y ^] 


d [pc] 


{B-V) 


{J-K) 


sp. type 


TY CrA 
cmpl 
cmp3 
cmp5 


7421-1126 
7421-1061 
7421-1163 


176423 
176497 


210829 
210834 
210840 


-0.2(1.3) 
6.1(1.4) 
3.3(1.3) 


-28.6(1.2) 
-27.3(1.4) 
-28.5(1.3) 


136^19 


0.463(32) 
0.291(20) 
0.170(10) 


0.813(23) 
0.240(23) 
0.080(28) 
1.939(33) 


B8e 

AlIII/IV 
AOIV 



Table 2. Overview of telescopes/instruments and detectors used to obtain photometry of TY CrA. Abbreviations of the observatories 
(Obs.): LS - La Silla (29°15f0 S, 70°44f0 W); CA - Cerro Armazones (24°36.'0 S, 70°11.'0 W); CT - Cerro Tololo (30°10.'0 S, 70°48f0 
W). 



Obs. 


Telescope /Instrument 


Aperture 

[cm] 


Filters 


Detector 


Size 

[pix] 


FoV 


LS 


NTT/SOFI 


358 


JHKs 


Hawaii HgCdTe 


IKxlK 


4f66 X 4f66 




REM/REMIR 


60 


JHKs 


Hawaii I 


0.5KX0.5K 


lO.'O X lO.'O 




REM/ROSS 


60 


VRI 


Apogee Alta 


IKxlK 


lO.'O X lO.'O 


CA 


VYSOS6 


15 


BVRI 


Apogee Alta U16M 


4Kx4K 


2°42.'0 X 2°42f0 


CT 


1.3m/ANDICAM 


130 


BVRI 


Fairchild 447 


2Kx2K 


6.'0 X 6f0 






130 


YJHK 


Hawaii HgCdTc 


IKxlK 


2.'4 X 2'A 



-0.1 




0.9 I 1 , 1 , 1 , 1 , 1 , 1 , L_ 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 

Phase 



Figure 1. The JH LC of TY CrA obtained by us with REMIR 
(2009 and 2011) and SOFI (2006) at La Silla. The LC is incom- 
plete because of bad weather and pointing problem. The primary 
minimum was obtained with SOFI in 2006 and the rest of LC 
was observed with REMIR in 2009 and 2011. The observations 
from 2011 were reduced with respect to HD 176423. For other ob- 
servations we used cmp5 as the comparison star (see Table [TJ. 
Therefore, to obtain whole LC we took into account the bright- 
ness difference of the comparison stars (see the text). Note that 
the errors of LC points in JH photometry are within the range 
of 0.003 - 0.010 mag. 



The observations during the first night suffered from 
strongly variable conditions and were not used in the anal- 
ysis. The primary minimum observed in the second night is 
of good quality. In addition to precise minima timing (see 
Section 3), the data were used for the LC analysis. All SOFI 
observations use cmp5 as the comparison star (see Table [1] 
and Figure [2]| 



Additional infrared photometry w as taken in 2009- 
2011 with RE I VQ (Rapid Eye Mount, IZerbi et all I2OOII : 
lAntonelli erall 1 20031 : IZerbi et all |2004| ) (60cm diameter 
rapid response telescope) operated at La Silla, Chile. The 
telescope hosts two instruments: REMIR, an infrared imag- 
ing camera, and ROSS, a visible light imager. The two cam- 
eras can observe simultaneously thanks to a dichroic plate 
placed before the telescope focus. 

In order to fully cover the NIR JHR^ LC we got REM 
observing time by three proposals (in the periods A0T19, 
AOT21 and AOT23). The observations strongly suffered 
from bad weather and technical problems. The observing 
programme has been completed to a low percentage even 
though additional DDT time was granted. The first obser- 
vations taken in AOT19 and 21 followed the standard ap- 
proach of differential photometry, i.e. a pre-defined field is 
observed containing the target and the selected compari- 
son star (cmp5) (see Figure [21). However, the REM frames 
are affected by pointing and off-centering problems, so that 
comparison star and sometimes also the target were off the 
frames. The 2010 data (AOT21) are of low quality and were 
excluded from the LC analysis. For period AOT23, we fol- 
lowed a different observing strategy. Since the pointing was 
not sufficiently precise, we decided to ob serve other fields in - 
clud ing compa ri son st ars already used bv lCasey et al.l (|l998l ) 
and IVaz et al.l (|l998l ). and to center the pointing at these 
stars (HD 176423 (cmpl) and HD 176497 (cmp3)). These 
fields and TY CrA were observed in an alternating sequence 
with high frequency which must not be lower than the fre- 
quency of variations of the near-infrared background. The 
observational strategy was further improved by observing 
selected critical phases (ingress, egress, minimum) of many 

^ http://www.rem.inaf.it/ 
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eclipses distributed over different nights instead of observing 
only a few eclipses continuously from ingress to egress. This 
way, the loss of a single night due to weather or technical 
problems had less effect on the overall success of the obser- 
vations. Eventually, we obtained 2 secondary minima and 
out-of-eclipse observations around phases 0.25 and 0.75. 

The same observing strategy was used also in ad- 
ditional near-IR photometry carried out by ANDICAM. 
This is a simultane ous optical and NIR imaging instrument 
l|DePov et al.ll2003l) operated by the SMARTS consortium 
at the 1.3m CTIO telescope, Chile. The field of view of 
the NIR camera is roughly centered on the larger field of 
view of the optical arm. We observed in two nights in May 
2011; however, because of bad observing conditions (we ob- 
tained only several usable datapoints), this observation was 
excluded from the following analysis. 

Concerning comparison stars, note that, in the case 
of SOFI data (2006) and REM data from 2009, we used 
2MASS 19013729-3649210 (RA2000: 19:01:37.292, DEC2000:- 
36:49:21.04) as additional comparison star. The REM obser- 
vations in 2011 were measured with respect to HD 176423. 
The additional standard is fainter than HD 176423 in the 
J and H bands, by 1.799 ± 0.006 and 0.642 ± 0.004 re- 
spectively. The Ks band observations, especially in the case 
of REM, were very scattered. Hence, for the photometric 
analysis we could use only the JH LCs. 

The near-IR JH data reduction was implemented using 
IDlB, ESQ MIDA^, and ESQ ECLIPSeQ. The fiat-field 
frames and the bad pixel masks were created using the ESO 
ECLIPSE flat and average recipes. The fiat-field division, 
bad pixel correction, sky subtraction, and shift and add of 
the science frames was performed using the ESO ECLIPSE 
jitter recipel3 

In further reduction tasks, the IRAF packag^fjhas been 
used. In the next step the astrometric solution of the frames 
was performed. The WCS system was determined using the 
2MASS on-line catalogue for referenc^H. Then the aperture 
photometry of the target and comparison stars has been 
done. The instrumental magnitudes for several apertures ap- 
propriate for the seeing conditions at site were determined. 
The details of all near-IR observations are mentioned in the 
journal of observations in the Table O The differential ex- 
tinction was neglected. 

Note that the errors of the LC data points in JH pho- 
tometry are between 0.003 and 0.010 mag. In contrast to the 
optical images, there is no contamination due to the refiec- 
tion nebula. Figure [2] shows that in the near-IR there is no 
background emission from the nebula, contrarily to the op- 
tical where a strong contribution from the refiection nebula 
is observed. 




Figure 2. The field around TY CrA in the H passband (the image 
taken on June 27, 2009). 




Figure 3. The BVRI LC obtained by us with VYSOS6 at the 
Cerro Armazones in 2009. The errors of individual points in the 
BVRI LCs are in the range of 0.006-0.009 mag. The LCs were 
arbitrarily shifted in magnitudes for clarity. 



^ http:/ /idlastro. gsfc.nasa.gov/contents. html 
^ http:/ /www. eso.org/sci/software/esomidas/ 
^ http:/ /www. eso.org/sci/software/eclipse/ 

^ Each REMIR jitter sequence consists of five frames. The jitter 
pattern has a pentagon shape and is exactly the same for each 
sequence but might be rotated. Mostly, the offsets were made 
available to the recipe beforehand and refined by cross-correlation 
search when possible. In many cases, however, the jitter recipe 
even succeeded in performing a blind offset search. 
^ http://iraf.noao.edu/ 

http:/ /www. ipac.caltech.edu/2mass/releases/allsky/ 



Interestingly, we noticed an asymmetry in the primary 
minimu m in the near-IR S OFI data which was already no- 
ticed bv lCasev et al.l (|l998l ) in the optical. 

2.2 5 ra/ CCD photometry 

The optical CCD photometry of TY CrA was performed by 
VYSOS6, a wide field imaging instrument locate d at Cerro 
Armazones Observatory, Chile l|Haas et al.ll20lil l. 

Observations of TYCrA were taken in 2009 and 2011 
(Fig. [Sjl. By the second campaign in 2011, VYSOS6 got an 
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Figure 4. Primary minimum in BVRI obtained with VYSOS6 
at the Cerro Armazones observatory in 2011. 

identical twin installed on the same mount. One telescope of- 
fers the filter combination VR and the other BI , so that ob- 
servations are taken in two filters simultaneously. Both tele- 
scopes are equipped with an Apogee Alta U16M 4096x4096 
pixels CCD camera each, providing a 2°42'.0 x 2°42'.0 field 
of view. During the first observing campaign we have cov- 
ered the whole LC in BVRI, used for the analysis described 
in the following. Unfortunately, in the case of the 2011 pho- 
tometry, we obtained only an incomplete primary minimum 
(see Fig. [Jjl . The secondary minimum was not observed at 
all because of bad weather. 

For the optical diff'erential photometry, we also used 
HD 176423 and HD 176497 as comparison and check star, 
respectively. The stars were found to be constant within 
the observational accuracy throughout the observing period. 
We reduced all optical data in a standard manner (bias, 
dark, fiat-field corrections), and aperture photometry was 
performed using IRAF and custom written tools. 

The LCs obtained in 2009 with VYS0S6 show that the 
depth of the secondary minimum (at phase 0.5) increases 
towards the redder bands {R, T) and reaches almost 0.10 mag 
in the /-band. The errors of individual points in the BVRI 
LCs are in the range of 0.006-0.009 mag. 

The visual LCs of TY CrA are contaminated by the re- 
flection nebula. Part of its light is subtracted as background 
within the aperture photometry. Any remaining contribu- 
tion is corrected for by fitting the third light in the LC anal- 
ysis. 



Table 3. The journal of CCD (BVRIJH) observations of TY 
CrA obtained at La Silla (REMIR, SOFI) and Cerro Armazones 
(VYSOS6) observatories in the period 2006 - 2011. The Julian 
date HJDmean is corresponding to the beginning of observation. 
Nim. — number of CCD frames in R and J filter. 



Date 




Phase 


Filter 






2 400 000+ 










VYSOS6 


Jun 26, 09 


55009.719 


0.184 - 


0.211 


BVRI 


24 


Tun 97 no 
J un z / , uy 




n A7fi 
u.^ / u 


u.oyo 


j3 V m 




Jun 28, 09 


55011.786 


0.899 - 


0.933 


BVRI 


25 


Jul 22, 09 


55035.644 


0.158 - 


0.181 


BVRI 


8 


Jul 23, 09 


55036.505 


0.456 - 


0.565 


BVRI 


98 


Jul 24, 09 


55037.648 


0.852 - 


0.882 


BVRI 


20 


Aug 26, 09 


55070.530 


0.234 - 


0.263 


BVRI 


23 


Aug 27, 09 


55071.518 


0.576 - 


0.609 


BVRI 


17 


Aug 28, 09 


55072.492 


0.914 - 


0.993 


BVRI 


44 


Aug 31, 09 


55075.494 


0.953 - 


0.031 


BVRI 


59 


Sep 01, 09 


55076.559 


0.321 - 


0.349 


BVRI 


25 


Sep 02, 09 


55078.496 


0.992 - 


0.048 


BVRI 


41 


May 27, 11 


55709.812 


0.699 - 


0.711 


BVRI 


20 


May 29, 11 


55711.851 


0.406 - 


0.420 


BVRI 


20 


May 30, 11 


55712.674 


0.690 - 


0.703 


BVRI 


19 


Jun 04, 11 


55717.879 


0.492 - 


0.506 


BVRI 


20 


Jun 06, 11 


55719.707 


0.125 - 


0.172 


BVRI 


27 


Jun 08, 11 


55721.826 


0.858 - 


0.874 


BVRI 


21 


Jun 09, 11 


55722.674 


0.152 - 


0.209 


BVRI 


40 


Jun 10, 11 


55723.681 


0.501 - 


0.514 


BVRI 


21 


Jun 11, 11 


55724.809 


0.891 - 


0.906 


BVRI 


20 


Jun 12, 11 


55725.893 


0.191 - 


0.266 


BVRI 


29 


REM-IR 


Apr 02, 09 


54923.907 


0.478 - 


0.481 


JH 


5 


Jun 13, 09 


54995.740 


0.344 - 


0.368 


JH 


29 


Jul 23, 09 


55036.469 


0.443 - 


0.456 


JH 


12 


Jul 24, 09 


55036.509 


0.457 - 


0.489 


JH 


39 


Jul 25, 09 


55037.801 


0.344 - 


0.368 


JH 


29 


Aug 20, 09 


55063.585 


0.830 - 


0.852 


JH 


28 


Sep 05, 09 


55079.567 


0.363 - 


0.384 


JH 


24 


Jun 12, 09 


55359.749 


0.352 - 


0.380 


JH 


24 


Jun 22, 09 


55369.732 


0.808 - 


0.832 


JH 


21 


Apr 25, 11 


55676.713 


0.074 - 


0.126 


JH 


9 


May 01, 11 


55682.712 


0.150 - 


0.494 


JH 


6 


Jun 23, 11 


55735.680 


0.486 - 


0.539 


JH 


12 


Aug 14, 11 


55787.552 


0.442 - 


0.527 


JH 


4 


SOFI 


Sep 07, 06 


53986.460 


0.964 - 


0.042 


JH 


242 



3 MINIMA DETERMINATION AND THE 
LIGHT-TIME EFFECT 

Our observations enabled us to determine 3 minima 
times of th e TY CrA eclipsing paiiQ, which add to the 
Casev et al. (1998') data and one recent observation by 
Paschke (2010). They are listed in Tab. H The times of 
minima were determined sep arately for all filters using the 
iKwee fc van WoerdenI l|l956l ) method and then the weighted 
averages were calculated. 

^ Partially covered eclipses were not used for the timing. 



For computation of minima times we have used only 
data in the phase interval ±0.05 around either of the min- 
ima. This ap proach minimizes the infiuence of the minima 
asymmetries ijvan 't Veerll 19731 ). 

Figure[5]shows the (O-C) residuals for all available min- 
ima (Table |4]) with respect to the optimum ephemeris de- 
termined by the LC modeling (Table [5]). As expected, for 
a detached binary, the orbital period of the system is sta- 
ble. Nevertheless, the (O-C) residuals exhibi t deviations sig- 
nifican tly larger than their standard errors. ICorporon et al.l 
(| 19961 ) presented five different spectroscopic orbital solutions 
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Table 4. Primary and secondary minima times of TY CrA. For 
some mi nima, the filters us ed were not given in the original pub- 
lication jCasev et a l.' 1998). For new minima observed in more 
than one filter the weighted average is given. 



HJD 




Filter 


Type 


2 400 000+ 


[days] 






Casev et al. C 1998) 


47694.7971 


0.0009 




pri 


47700.5740 


0.0005 




pri 


47710.68 


0.08 




sec 


48783.8703 


0.0013 




pri 


49153.6302 


0.0032 


y 


pri 


49160.85 


0.02 


y 


sec 


49163.74 


0.01 


yb 


sec 


49514.7310 


0.0003 




pri 


49527.738 


0.005 




sec 


49537.8414 


0.0008 




pri 


49543.6193 


0.0003 




pri 


This work 


53986.55284 


0.00038 


JH 


pri 


55036.64046 


0.00035 


IR 


sec 


55075.63231 


0.00016 


BVRI 


pri 


Paschke ("2010) 


55387.619 


0.005 


V 


pri 



for periods of the outer orbit P3 — 126-270 days. The most 
probable are their solutions #1 (P3 = 270 days) and #2 
(P3 — 268 days). The expected light-time effect (hereafter 
LITE) amplitude modulating the binary eclipse timings was 
computed for solution #1 and the following masses of the 
components: Mi = 3 Mq, M2 = 1.6 Mq and M3 = 1.2 Mq. 
Using the published major axis a = 1.47 a.u., the outer or- 
bit inclination i — 20°, and the above masses we get the 
projected semi-major axis of the eclipsing pair as: ai2 sin i 
= 0.104 a.u. corresponding to 51.9 light seconds (or 0.0006 
days). In the view of the most precise minima having a — 
0.0002 days and probable systematic errors (caused e.g. by 
activity of TY CrA), the detection and useful analysis of the 
timing signal caused by the third component can hardly be 
performed using the ground-based photometry. 

Because of the short period of the outer orbital pair the 
RV observations are more sensitive and useful in the decisive 
determinat ion of the orbital elements. To summarize, we 
agree with lCasev et all (Il998l ) in that the photometric O-C 
diagram of the inner eclipsing binary cannot be explained 
by orbital sol utions for the sp e ctros copic tertiary suggested 
previously bv lCorporon et al] l|l996[ ). 



4 LIGHT CURVE ANALYSIS 

Unfortunately, optical {BVRI) data and near-IR (JH) pho- 
tometry were obtained in different time intervals (see Ta- 
ble |3]). Most of the optical photometry was taken in 2009, 
while the near infrared data used for analysis in 2006, 2009, 
and 2011. Because of the extreme stability of LCs, we at- 
tempted to model all photometric data simultaneously. Sep- 
arate solutions of visual and NIR data were also performed. 
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* Caseyetal. 1998 
° This work 

• Paschke 2010 
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5 





43000 50000 52000 54000 56000 



HJD 2400 000 + 

Figure 5. (O-C) diagram for all published and new minima (ob- 
tained from optical and NIR LCs) timings corresponding to the 
optimal ephemeris obtained by the simultaneous multi-colour LC 
analysis (Table |5)l. 



Our photometry of TY CrA is the sum of the eclipse LC 
of the inner binary and a constant light contribution of the 
unresolved third component and the visual companion (0.3 
arcsec from the unresolved triple). Moreover, TY CrA and 
HD 176386 (1^1 southwest of TY CrA) illuminate the reflec- 
tion nebula NGC 6726/6727, whose light is a sign ificant con- 
tributor to the observed flux (jCasev et al.lll998h . The light 
from the two stars in the eclipsing pair is of principal inter- 
est here. The remaining sources produce a constant offset, 
so-called third light, decreasing the LC amplitude and com- 
plicating the analysis. Because of the three different sources, 
the third light cannot be modeled as a single stellar source 
(deflned by its fractional radius, and effective temperature). 
Hence the third light will be treated without any assump- 
tions and independently in the individual passbands. 

The latest version of the code ROCHE was employed 
(|Pribullall20il ). The LC of the system was modeled assum- 
ing (1) radiat ive envelopes w ith bolometric albedos Ai — A2 
= 1.00 (e.g. iRucinskH Il969f l and gravity darkening coeffi- 
cien ts q\ = Q2 = 1-00 appropria te for the B8V spectral type 
fsee IPantazis fc Niarcho3ll998l ). (2) limb darkening coeffi- 
cients automatical ly interpolated from extensive tables of 
Ivan Hammd l|l993l ) for average surface gravity logg, aver- 
age Te/ / and given passband, (3) solid body rotation of the 
components (no differential rotation), (4) synchronous rota- 
tion of the secondary component F2 = 1 and no rotation of 
the primary Fi = (see below) , (5) rotational axes perpen- 
dicular to the orbital plane, (6) circular orbit, (7) no pho- 
tospheric spot(s), and (8) local intensity computed using 
model atmospheres (local logg, Te//, and A) for the solar 
chemical composition. The mutual irradiation effects were 
approximated by single reflection. Surface grids were derived 
from a icosahedron in such a way that surface elements are 
close to equilateral triangles of similar size. 

In the optimization process the following parameters 
were adjusted: To,P - corresponding to the primary min- 
imum and orbital period, i - inclination angle, f2i,f22 - 
generalized surface equipotentials, T2 mean temperature of 
the less massive and colder secondary component, + 
hi^) - passband-specific unit-less normalization factors. 
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and l-i{X) - third ligh10. The temperat ure of the p r ima: 
Ti = 12000±500 K was adopted from ICasev et al.l jm 
while the spectroscopic orbit {Ki — 85.2 ± 0.2 km.s~^, 
K2 — 164.6 ± 1.6 k m.s~^) was adopted (not adjusted) from 
ICasev et al.l (|l995l ). 

Before arriving at the final solution we extensively 
tested the subsynchronous rotation hypothesis. Because of 
the small fractional radiu^3 of the primary component 
(r2 ~ 0.1) the asynchronism has a very small effect on the 
LC shape. The relative change of the equatorial radius of 
the primary component between the case of no rotation {Fi 
= 0) and synchronous rotation {Fi = 1) is about 0.8% only. 
Hence we assume that the primary does not rotate at all. 
The multi-colour LCs were solved by the method of differen- 
tial corrections. The optimisation process always resulted in 
the same parameter set (we used different starting parame- 
ters) indicating uniqueness of the solution. 

The derived photometric elements are presented in Ta- 
ble El The values found for I3 are 0.8965, 0.6969, 0.6791, 
0.6940, 0.3032, and 0.4205 for B, V, R, I, J and H, re- 
spectively, showing the increasing influence of the nebular 
contribution for shorter wavelength (see Table [5}, the addi- 
tional stellar component(s) dominating in the NIR. The best 
fit to the BVRIJH LCs is displayed in Fig. |6l In the case of 
the //-band data the primary minimum together with fit 
has lower {^^ 0.08 mag) amplitude than in the other filters. 
This fact is due to the bad data quality in the primary min- 
imum. The scattered points have been removed so that the 
amplitude of the fit is lower and the third light contribution 
is larger. The BVRIJH LCs solution for synchronous rota- 
tion and rotational axes perpendicular to the orbital plane 
predicts vi sini = 25.4 km s~^ and V2 sini — 40.0 km s~^. 

The absolute parameters of the primary and secondary 
components of TY CrA (corresponding to the simultane- 
ous fit to the (BVRIJH) LCs) are given in Table El The 
values of most parameters, ex cept of Ri and Li , are con- 
sistent with those derived by ICasev et all (|l998l ). We ob- 
tained about 19% smaller diameter of primary and about 
10% larger diamet er of the seco ndary component compared 
to values given bv lCasev et al. 1 11998) . Such a difference in 
diameter translates to even larger difference in luminosity: 
our value for the primary is by 42% lower. This discrepancy 
will be discussed in the following section. 

T o better compare our results and those of IC asev et al.l 
l|l998l) , we also performed an analysis of the optical BVRI 
LCs only. The resulting parameters giv en in Table [6] are sti ll 
not consistent with those presented bv lCasev et al.l l|l998l ): 
Ri and Li are about 15% and 35% smaller, respectively. 

In addition, we solved the LCs in the BVRI passbands 
without t hird light reproduc ing very closely the parameters 
given by ICasev et al.l (| 19981 ). However a good solution in 
the J and H bands cannot be achieved without admitting a 
third light contribution. The conclusion is that the near-IR 
data allows us to constrain the third light very well. 



"T ' 1 ' 1 ' r 




^ The third light is expressed in ROCHE as = L-^/^L^ 
relative to the major axis a, unitless quantity 



-L2). 



Figure 6. Best fits to the observed BVRIJH LCs of TY CrA 
resulting from the simultaneous modeling of all bands. The orbital 
phases are determined according to the new ephemeris given in 
Tabled 



5 THE EVOLUTIONARY STATE OF TY CRA 
5.1 Mass-radius diagram 

Masses and radii are the most fundamental parameters de- 
rived from the observations. They are hardly affected by 
uncertainties in the adopted effective temperature of the pri- 
mary component. We varied the primary temperature in the 
light curve analysis from 10,000 K to 14,000 K in order to see 
whether there is any effect on the resulting masses and radii. 
While luminosity varies drastically from 30 to more than 
70Z/Q, the changes in masses and radii remain well below 
the formal error bars. In other words, the derived values of 
mass and radius are very robust w.r.t. an uncertain primary 
Tcff . 

The radii are indicative of the evolutionary state of 
both components. The mass-radius diagram in Fig. [Tlcom- 
par es the loci of both components to previous assessments 
by ICasev et al.1 (Il998h and precise data of other bina- 
ries llTorres et al.ll2O10D . In addition, PISA isochrones of 
iTogneU i et al. | (I2OIII) of differe nt age and metallicity are 
overplotted. ICasev et all (|l998l ) have already noticed that 
TY CrA has reached the main sequence while the secondary 
still resides in the pre-main sequence phase. 

Not only the secondary but also the primary is an im- 
portant object since it is located in an otherwise empty re- 
gion of the mass-radius diagram at small radii. Our value 
of the primary radius is even lower than the value given 
bvlCasev et al.l (| 19981 ). The comparison to the [Torres et al.l 
l|2010l ) data shows that the primary component seems con- 
siderably undersized. The primary only agrees with the 
models when adopting a metallicity of Z=0.003 or less 
(Figs. Oj and [71;). For comparison. Fig. (Tji shows the 
Z =0.0125 isochro n e acc o rding to the metallicity adopted 
by ICennaro et al.l l|2012l ) . ICasev et al.l l|l998l ) compared to 
models at even higher metallicity. Sun-like and Hyades-like. 
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Figure 7. Mass-radius diagram. The new measurements of 
T Y Cr A (c ircles) are compared to the previous results obtained 
bv lCasevet al. { 199^, squar es) and the data of other binaries com- 
piled bv Plbrres et aU 1 I2OIOI ). The measure ments are compared to 
theoretical isochrones for 3 and 5 Myrs from lTognelli et al.! boil l) 
at different metallicity: a) Z=0.0125, ^=0.255, Xd = 2 ■ 10"^, 
a = 1.20; b) Z=0.001, y=0.232, Xd = 4 • IQ-^, a = 1.20; c) 
Z=0.003, y=0.236, Xd = 4 ■ 10"^, a = 1.20. 



Table 5. Photometric elements (coming from BVTJ/and Ji? si- 
multaneous analysis) of the eclipsing pair in the TY CrA system 
and their standard errors, cr (number in parentheses gives the er- 
ror of the last digit) - i - inclination; e - eccentricity; q = m2/mi 
- mass ratio an d K\+K2 - sum o f radial- velocity semi-amplitudes 
(adopted from lCasev et a"l]|l995l ); f2i^2 - generalized surface po- 
tentials; /i^2 - asynchronous rotation factors, T2 - effective tem- 
perature of the secondary component; rp^ii '"sec - volume mean 
fractional radii; /{L^ -|-L^) - monochromatic lights and (3 - 
contribution of the third light. Parameters fixed (not adjusted) 
in the solution are denoted by a superscript " /" . The uncertainty 
of the secondary-component temperature contains quadratically 
added 200 K error resulting from the 500 K error of the primary- 
component temperature estimated by Ca sey et al.l l|l998l l . The er- 
ror has not been propagated to uncertainties of paasband-specific 
relative lights Li/(Li -|- L2)- 



Parameter 


TY CrA 


HJDo 


2 449 537.8423(25) 


P [days] 


2.88877912(21) 


i [°] 


85.1(10) 


e 


0^ 


Hi 


10.4(10) 


Q.2 


4.5(6) 


Fi 


o.o-f 


F2 


1.0^ 


1 


0.518^ 


Ki +K2 [km.s-l] 


249.8^ 


Ti [K] 


12000-'' 


T2 [K] 


4900(250) 




0.102(10) 




0.161(9) 


Lf/(Lf +Lf) 


0.9674(14) 


LY/{LY + LX) 
L^/{L^ + L^) 
L{/{L{ + Li) 


0.9177(32) 


0.8581(10) 


0.7905(9) 


L(/{L( + Li) 


0.6592(13) 


0.5400(11) 



'3 

'3 



0.8965 
0.6969 
0.6791 
0.6940 
0.3032 
0.4205 



The secondary agrees with a very young age of some 
3-5 Myrs for any set of isochrone s, in accordance with 
previous estirnates b y I Casey et al.l 
iGennaro eraTI l|2012l . « 4 - 5Myrs). 



3 Myrs) and 



5.2 HR diagram 

The evolution of the radiative quantities luminosity and ef- 
fective temperature is studied in the HR diagram (Fig. [S]). 
Again, the new measurements are comp ared to the previou s 
location in the HR diagram according to lCasev et al.l (|l998l ) 
and to PISA models, now including stellar trac ks. For clarity 
and si nce we do not gain additional insight, the I Torres et al.l 
(|201(j ) data is not shown in Fig. [H] 

Figure [S^ aims at reproducing fig. 8 of iGennaro et al.l 

who fitted the location in the HR diagram based on 
the Casey et al.l (|l998l ) data. However, Fig. [8^ looks more 
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Figure 8. HR diagram. The new measurements of TYCrA 
(circles) are com pared to the previous results obtained by 
ICasev et al] l ll998l. squares). The m easurements are compared to 
PISA models llTognelli et al]|201ll) at different metallicity. Theo- 
retical isochrones for 3 Myrs and 5 Myrs are displayed as well as 
mass tracks for 1.6, 1.7, 3.0, and 3.2 Mg: a) Z=0.0125, Y=0.255, 
Xd = 2 • 10-5(solid)/4 • 10"^ (dashed), a = 1.20 (sets with 
differe nt are used to enable comparison with lGennaro et al.l 
l|2012l ). and to show that different Xd has relatively small effect); 
b) Z=0.001, y=0.232 (solid) / 0.251 (dashed), Xd = 4 ■ IQ-^, 
a = 1.20 (sets with different Y are used to enable comparison with 
other metallicities, and to show that different Y has relatively 
small effect); c) Z=0.003, ^=0.236, Xq = 4 ■ IQ-^, a = 1.20. 



Table 6. Absolute parameters of the TY CrA eclipsing pair. The 
second column corresponds to solution #1 from the analysis of 
BVRIJHhCs. The third column shows solution #2 resulting from 
analysis of BVRI LCs only. The mass ratio, q and Ki + K2 were 
adopted from the circular solution of lCasev et al.l lll99d). All re - 
sults are compared with values published bv lCasevet''al] lll998l) . 
The larger masses for solution 2 result from lower inclination an- 
gle, i = 82.9°. 



Par am. 


Solut. 1 


Solut. 2 


Casey et al. 




[BVRIJH] 


[BVRI] 


(1998) 


Mi/M© 


3.11(7) 


3.15(7) 


3.16(2) 


M2/MQ 


1.61(2) 


1.63(2) 


1.64(1) 


R1/R0 


1.46(15) 


1.53(13) 


1.80(10) 


R2/RQ 


2.30(14) 


2.21(12) 


2.08(14) 


Ti [K] 


12000(500) 


12000(500) 


12000(500) 


T2 [K] 


4894(250) 


4977(450) 


4900(400) 


Li/Lq 


40(10) 


44(10) 


67(12) 


L2/LQ 


2.7(7) 


2.7(10) 


2.4(8) 


log 51 [cgs] 


4.60(8) 


4.57(7) 


4.43(6) 


log S2 [cgs] 


3.92(5) 


3.96(4) 


4.02(5) 


a [Rq] 


14.29(9) 


14.35(9) 


14.29(9) 



complicated since the very same tracks and isochrones are 
not available to us. Therefore, the two mass tracks clos- 
est to the measured masses are shown for each compo- 
nent of TY CrA. In addition, models for different deuterium 
abundances are shown in order to be able to compare with 
Figs. (Sja and [8};. However, the effect of deuterium abun- 
dance is small compared to the variations due to age, mass, 
and metallicity. Obviously, the tracks and isochrones are at 
variance with the new primary data. An isochrone of older 
age might still be consistent with the primary, but would 
disagree with the secondary. Furthermore, the discrepancy 
with the primary in the mass-radius diagram (Fig. [7^) would 
not be improved. 

In contrast to the mass-radius diagram, the overall sit- 
uation in the HR diagram does not improve when choos- 
ing low-metallicity models. At least, choices of low metal- 
licity allow for agreement of the primary with isochrones of 
3 Myrs (Fig. 111313 and 5 Myrs (Fig. |8l;), respectively. The 
disagreement with the mass tracks, however, is exacerbated. 
Moreover, the low-metallicity tracks disagree with the sec- 
ondary in the HR diagram. For Z^O.003 (Fig. [8};), the evo- 
lutionary tracks agree at best marginally with the new sec- 
ondary parameters. Moreover, no isochrone can be found 
which matches both components. In summary, low metallic- 
ity does not explain the small primary radius, even though 
this is suggested by the mass-radius diagram. 

It is interesting to note, that the models of Z=0.0125 
might agree with the observations in the HR diagram when 
adopting a different value of the primary effective tempera- 
ture which is an input parameter in the LC analysis. While 
there is no significant effect on the radius measurements, 
the luminosities derived vary substantially. Nevertheless, the 
discrepancies in the mass-radius diagram will persist. 

In Fig. [S]], models for two different helium abundances are 
shown to allow for a better comparison with Figs. |8^ and [8};. 
Again, the effect is comparably small, in a way similar to the 
effect of varying deuterium abundance. 
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6 SUMMARY 

We present new optical and NIR photometry of TYCrA 
taken in 2006-2011. The BVRI photometry t aken comple- 
ments the Stromgren photometry studied by ICasev et al.l 
l|l998l l. For the first time, we present NIR light curves of the 
secondary minimum taken in the J and H bands. In con- 
trast to the optical light curves the secondary minimum is 
much better revealed in the NIR and is more than 0.1 mag 
deep 

We compiled all minima times available to us so that 
the O-C diagram of TYCrA now spans over two decades. 
The light-time effect of the spectroscopic tertiary is too 
small to explain the changes in the O-C diagram. There are 
probably other effects complicating the analysis like shifts 
of the observed minima due to photometric out-of-eclipse 
variability. We find, that RV measurements are more suit- 
able to solve the tertiary orbit. One requires high-resolution 
high S/N spectroscopy with an echelle spectrograph at a 
2-m class telescope covering at least one orbital period of 
the unresolved triple system (2 70 days most p robably) . The 
Fourier-domain technique (see iHadraval Il995l ) then can be 
used to disentangle the component's spectra for further anal- 
ysis (Toff, logg, log [m/X]) and to improve the spectroscopic 
orbits. The orientation of the components' rotational axes 
could be found analyzing the McLaughlin rotational effect 
during eclipses of the inner pair. Long-baseline interferome- 
try would help to determine orbital parameters of the outer 
orbit of the triple and to determine the distance to the 
system (comparing interferometric and spectroscopic semi- 
major axis of the inner triple system) and the total mass of 
the triple. 

For the first time, we present a multi-band solution of 
the light curve of the inner eclipsing binary in both opti- 
cal and neax-IR bands. The near-/i? data requires to ac- 
count for a third light contribution and thus also allows us 
to constrain the colours of the third light. The third light 
is composed of at least three sources which are the reflec- 
tion nebula, the spectroscopic tertiary, and the visual fourth 
component. 

Only the radius and the luminosit y of the primary d o 
not agree with the values published bv lCasev et al.l l|l998h . 
Agreement is not improved by restricting the analysis to the 
optical data. While the primary luminosities could be recon- 
ciled by adjusting the effective temperature of the primary, 
the derived primary radius is very robust and would not 
change. 

Indeed, the absolute parameters of TY CrA do not agree 
with expectations. The primary component seems under- 
sized when compared to other stars of same mass and to 
recent PISA evolutionary tracks and isochrones. The HR di- 
agram shows that low metallicity cann ot explain the sma ll 
radius of the primary. As speculated bv lCasev et al.l (|l998h . 
the subsynchronous rotation of the primary and effects on 
its internal structure might play a role. 

A very interesting question is whether the fourth com- 
ponent is actually bound to the TY CrA system. In order to 
decide whether the fourth component is bound, there is need 
for single-epoch AO observation of the system to establish 
the physical bond of the fourth visual component O'.'S away 
in 2002. If it was a background star, the proper motion of 
TY CrA (see Table [1]) would make the visual companion 



unresolved around 2012. Because of the high total mass of 
the whole system the orbital period of the outermost visual 
orbit is probably rather short: assuming a distance of d = 
130 pc, the 0'.'294 separation observed in 2002 as the semi- 
major axis and the total mass of the system as 7 Mq, we 
get Pi234 = 89.3 years. 
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